A Dilute Ising Ferromagnet on a Hierarchical Lattice with Attractive Biquadratic 

Interactions 
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This paper considers a dilute Ising ferromangnet with annealed vacancies and attractive bi- 
quadratic interactions. Phase diagrams have been calculated while varying the temperature and 
concentration of annealed vacancies in the system while maintaining constant, attractive biquadratic 
couplings. These results have been produced using renormalization group analysis with a hierarchi- 
cal lattice. Critical exponents have been calculated and each basin of attraction interpreted. 
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I. INTRODUCTION 

The Blume-Emery-Grittiths (BEG) model [H is a 
spin-1 Ising model that is useful in probing the nature 
of systems with fluctuation in magnetization and den- 
sity. The Hamiltonian, shown in Equation 1, has the 
bilinear {Jij), biquadratic {Kij) and crystal- field inter- 
action (Aij) terms as shown in the Hamiltonian in Eq. 



j3H 



„2„2 



5^ JijSiSj + Ki- 

with Si = 0, ±1 



(1) 



With these three coupling coefficients we can probe 
ordering in various systems as we vary temperature (~ 
1/J), concentration (~ A/J) of annealed vacancies, and 
the clustering bias (~ K/ J) . 



J2l,,{s^s,+s,s^^) (2) 



In addition to the bilinear, biquadratic and crystal- 
field interaction terms we must also consider odd sector 
contributions to the Hamiltonian, as shown in Eq. (2). 
These contributions must be included in order to obtain 
a complete description of higher order phase transitions 
arising in our system. 

It should be noted that each summation in Eqs. 1 and 
2 is over nearest- neighbor (ij) pairs of our lattice unit 
structure including the magnetic (H) and crystal-field 
(A) interaction terms. The summation for each of these 
terms has been shifted from sites to bonds for computa- 
tional expediency with the overall effect being that A in 
Eq. (1), and H in Eq. (2), is the chemical potential, and 
magnetic field, per bond divided by two. 

The introduction of density as an added degree of free- 
dom in Ising systems was originally proposed to investi- 
gate the superfluid transition in He^ — He^ mixtures [l[ . 
Since, these spin-1 Ising systems have been extended and 
used to explore structural glasses 0, binary alloys, mi- 
croemulsions [3|, materials with mobile defects, binary 




fluids, aerogels Q, frustrated percolation systems 
among many others. 

Fixed points, critical phenomenon and resulting phase 
diagrams can be drastically altered due to underly- 
ing competing microscopic interactk)ns_in_yarious Ising 
systems. Mean-field methods 
renormalization-group techniques [lOI . 1111 . 11^ 
[Tg!] have been used in many previous studies to consider 
a range of qualitatively unique competing interactions 
using the Blume-Emery-Griffiths model. 

Renormalization-group methods have been employed, 
in conjunction with hierachical lattices, to investigate the 
effects of competing biquadratic interactions in a dilute 
Ising ferromagn et jl4l| , competing bilinear interactions in 
a BEG system [iS^, competing bilinear interactions [Tlj 
in a spin-1/2 Ising model, and simultaneous competition 
between crystal-field and biquadratic interactions in a 
BEG ferromagnet [T2|. 

Other renormalization-group probes yield tricritical 
and critical end point topologies linking first and second 
order phase boundaries for the case with K > Q. For neg- 
ative biquadratic coupling [K < 0), mean- field calcula- 
tions revealed two novel phases: one a high-entropy ferri- 
magnetic phase and the other displaying antiquadrupolar 
order, see reference 

The effects of attractive and repulsive biquadratic in- 
teractions on the global phase space was considered by 
Sellitto et al. using the replica symmetric mean- 

field approximation with quenched disorder in the bilin- 
ear interactions. A spin-glass phase was found with both 
first and second order transitions from the paramagnetic 
phase: the order of the transition largely dependent upon 
the crystal-field interaction. This study also found, for 
strong repulsive (K < 0) biquadratic interactions, an an- 
tiquadrupolar phase and at lower temperatures, an anti- 
quadrupolar spin-glass phase. 

Kabakcioglu et al. considered the effects of quenched 
random fields [13] and Falicov et al. the effects of 
quenched random bonds [l^ upon ordering, critical- 
ity and resulting phase diagrams. Random crystal fields 
were the focus of Branco et. al using real-space RG 
[Tsl . and mean-field approximations for both Blume- 
Capel and Blume-Emery-Griffiths model Hamiltonians, 
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respectively. 

The current study complements these earlier works as 
it considers a dilute Ising ferromangnet with annealed 
vacancies and attractive biquadratic interactions. The 
effects of temperature and concentration of annealed va- 
cancies upon ordering for a series of constant, attractive 
biquadratic couplings have been probed via calculation of 
several phase diagrams. Each phase diagram produced 
using renormalization group analysis with a hierarchical 
lattice. Critical exponents have been calculated for alU 
higher-order transitions and each basin of attraction in- 
terpreted. 




FIG. 1: The construction of an infinite hierarchical lattice by 
repeatedly replacing each nearest-neighbor interaction by the 
basic unit itself (Berker and Ostlund p^). 



II. HIERARCHICAL LATTICES AND 
RENORMALIZATION GROUP 

The basic recipe for generating an infinite hierarchical 
lattice from its basic unit is to repeatedly replace each 
nearest-neighbor interaction by the basic unit itself. The 
construction of a generic hierarchical lattice is shown in 
Figure 1. Figure 2 illustrates the construction of the 
hierarchical lattice [Hi [13] used for this study, with the 
difference being the complexity of the basic generating 
unit. 

Since hierarchical lattices yield exact renormalization 
group recursion relations for the coupling coefficients, 
phase diagrams can be produced and critical exponents 
determined very accurately. Therefore the results pre- 
sented in this study may be considered exact on the 
rather specialized lattice, or, they may be used as ap- 
proximations into these systems modeled on more re- 
alistic lattices. A range of very difficult problems has 
been subjected to study using hierarchical lattices as the 
medium. For example, random-bond 12 ill , random-field 
[ill, spin glass [ll,|2l|, frustrated [li!,ll3, [Ti] , directed- 
path H^l and dynamic scaling [2^ systems have all been 
investigated and better understood using these special- 
ized lattices. 

The renormalization-group solution for a hierarchical 
model, such as Figures 1 and 2, essentially reverses the 
construction process. Internal degrees of freedom are 
eliminated with each renormalization by summing over 
all configurations of the internal sites (represented by 
solid black dots in Figures 2a and 2b, and by ai,aj in 
Equation 6). 

With each renormalization, or rescaling, in our sys- 
tem we demand the partition function remain unchanged. 
This allows us to derive recursion relations that relate 
the exchange interactions at the two length scales. The 
new effective interactions J', K' , and A' are separated 
by a distance I' which is b lattice constants in the origi- 
nal system, where b is the length rescaling factor of the 
renormalization-group transformation. 




FIG. 2: Construction of the hierarchical lattice used in this 
investigation. (Reprinted from Journal of Magnetism and 
Magnetic Materials, 314, 69-74 (2007) D. P. Snowman, with 
permission from Elsevier) 

Ci,{J',K\A') = OiJ,K,A) (3) 
with /' = bl (4) 

Ci = Yl exp[-/3H] = J2 (5) 
with Ri{si, Sj) — ^ exp[— /3i7] (6) 

= 5]exp[-/3H'] = 5]i?K.s:,s;.) (7) 

with i?i(sj, Sj) — ^ exp[j'siSj 

+K'sy^-A'{sj + s'^) + &] (8) 
where G' is a constant used to calculate the free energy. 
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Equating the contributions to the partition function, 
Ri{si, Sj) and Ri {si, Sj) - corresponding to the same fixed 
configuration of end spins, Si,Sj and s^,Sj - at the two 
different length scales, I and I, allows us to calculate the 
actual renormalzation-group transformations. Thus re- 
lating the interaction strengths at the two length scales, 
I and I, can be derived: J'{J, K, A), K' {J, K, A), and 
A'{J, K, A). The derivation of these relations is pre- 
sented in more detail in Section 5. 

Parameter space is explored, phase diagrams mapped 
and transitions characterized using these recursion re- 
lations in conjunction with the initial values of J, K 
and A. Repeated iteration of the recursion relations 
from an initial starting point carries the system along 
a renormalization-group trajectories and eventually to a 
sink. 



III. PHASE TRANSITION 
CHARACTERIZATION 

Numerically differentiating the free energy density 
with respect to the appropriate variables yields values 
for the magnetizations, densities and nearest neighbor 
correlations for our system. The free energy density (di- 
mensionless free energy per bond ), /, can be expressed 
as 



/ ^ -— ^ ^6-"''G"(")(j("~i), if A("-i) (15) 



J' = Rj{J,K,A) (9) 
K' = Rk{J,K,A) (10) 
A'=R^{J,K,A) (11) 



Two types of components (see Fig. 2a and b) have been 
included in the basic unit, similar to pll. IT^. Iisl. 111. l26l| . 

for the calculations presented in this paper. One type 
of component, type A, allows internal spins to interact 
via a cross-link feature. The strength of the cross-link 
interaction is varied via the parameter p, with the case 
oi p — Q yielding the hierarchical model equivalent [l^ 
to the Migdal-Kadanoff [23,[11] decimation-bond moving 
scheme in two dimensions. 

A second type of component, type B, allows the end 
spins to interact via two different length connecting 
paths. As shown in Figure 2b, the first path has mi 
pairs of spins and the second m2 pairs of spins. The 
overall connectivity of the system is varied using two pa- 
rameters, PA and pb- These parameters, represent the 
number of unit structures, either type A or type B (see 
Fig. 2), used to form the basic generating unit for the 
hierarchical lattice as shown in Fig. 2c. The connectivity 
parameters, {p,mi,m2,PA,PB) — (4,8,9,40,1), used in 
the present investigation parallel those used in previous 
studies, see references [nl, [11 [H, H [13 , probing the 
effects of various types of competing interactions on a 
similar lattice. 



where F is the Helmhotz free energy and Nj, denotes 
the total number of bonds in the system. The free en- 
ergy density consists of a sum, over all iterations of the 
renormalization-group transformation, of the contribu- 
tions G"(") to the free energy density due to the degrees 
of freedom removed during each transformation. With 
each rescaling, or renormalization, the length scale of the 
system is reduced by a factor of b and the number of spins 
by a factor of b'^. 

The magnetization, m = ^ = ^j^, can be cal- 
culated by numerically measuring the shift in the free 
energy density with a small shift in the magnetic field, 
where Ng is the number of sites. Using a similar approach 
yields the density, p = . Nearest neighbor corre- 

lations of the bilinear, {siSj) — ^^'^ biquadratic, 

(sfs'j) — exchange interactions also aid in inter- 

preting resulting phases and transitions. 

Using the four thermodynamic quantities discussed 
above, transitions between the various basins of attrac- 
tion are characterized. First order transitions reveal 
themselves by discontinuities in densities, magnetiza- 
tions, or other first derivatives of the free energy. Critical 
transitions exhibit no such discontinuities in any of the 
order parameters. Exact expressions for the recursion 
relations allow us to calculate critical scaling exponents 
for our system. The next section explores this in greater 
detail. 



J* = Rj{J\K*,A*) (12) 



IV. RECURSION RELATIONS 



j{* — Rj^(^J* ^ K* , A*) (13) Fo'" each fixed end spin configuration, equating the 

contributions to the partition function, Ri(si,Sj) and 
Ri(si,Sj), from the two length scales - as introduced in 

^ ^ ^a{J ,K,A) (14) section 2 - allows us to write the following equalities for 

the type A structure shown in Figure 2a. 
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exp[-4A] + 2 exp[-2 J + 2K- A(6 + p)] 

+2 exp[2 J + 2K- A(6 + p)] + cxp[J(-4 +p)+ K{4 + p) - A(8 + 2p)] 

+2 cxp[- Jp + K(4 + p) - A(8 + 2p)] + exp[J(4 + p) + K{4 +p)- A(8 + 2p)] 

= exp[J' + K' - 2 A' + G] = Ri,[l, 1], 



(16) 



i?;[l,0] = exp[-2A]+2cxp[-J + if- A(4+p)]+2exp[J + ii:- A(4 + p)] 

+ exp[J(-2 +p)+ K{2 + p) - A(6 + 2p)] + 2 cxp[-Jp + K{2 + p) - A(6 + 2p)] 
+ exp[J(2 + p) + K{2 + p) - A(6 + 2p)] = exp[- A' + G] = Rv [1,0], 



(17) 



i?;[l,-l] = exp[-4A] + 4exp[2ii:- A(6+p)] 

+2 exp[- Jp + K{A + p) - Delta{8 + 2p)] 
= exp[- J' + if' - 2A' + G] = i?;. [1, -1], 



2 exp[ Jp + if (4 + p) - Delta{8 + 2p)] 



(18) 



ii;[0,0] = l + 4exp[-A(2+p)]+2exp[-Jp + ifp- A(4 + 2p)] 
+2 exp[Jp + Kp- A(4 + 2p)] = exp[G] = Ri^ [0, 0], 

I 



(19) 



Eqs. 16-19 can now be manipulated algebraically to 

yield expressions that relate the coupling coefficients be- 
tween the two length scales for the type A unit structure. 



Critical exponents can be extracted from the exact re- 
cursion relations above via linearization in the vicinity of 
the second-order transition under investigation. That is, 



, 1 i?r(l,l) 

^- = 2^°^ERi-rT) 



log 



Ri,{i,i)Ri'{i,-i)RUo,o) 



2 



i?f, (1,0) 
i?('(0,0) 



i?r(l,0) 
G;4 = logi?i.(0,0) 



(20) 

(21) 

(22) 
(23) 



The type B unit structure has recursion relations with 
the same form as in Eqs. 20-23, but the expressions 
(Eqs. 16-19) for the imdcrlying Rl{si,Sj) differ signifi- 
cantly. Combining the recursion relationships for both 
types of structures (type A and type B as shown in Fig. 
2), the renormalization relationships become 



K -K* 



Tjj{J-r)+TjK{K-K*) 
+Tja{A - A*), (27) 



Tkj{J-J*)+Tkk{K-K*) 
+Tka{A - a*), (28) 



A - A* = TAjiJ-J*)+TAK{K-K*) 

+Taa(A - a*), (29) 

where Tjj = Tkj = etc. and are evaluated 
at the fixed point in question. The critical relations above 
can be represented as a recursion matrix, with elements 
TxY and eigenvalues of the form 



J =paJa+PbJb 
k' = paK'j^+pbK'b 

A' =paA'^+PbAb 



(24) 
(25) 
(26) 



by 



where I = 2,4,6 



(30) 



where b is the length rcscaling factor (in our case b = 2) 
and yi represents the corresponding critical exponent for 
the l^^ eigenvalue. A similar linearization for the odd 
sector contributions {H and L) yields critical scaling ex- 
ponents 2/1 and 2/3. 
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Phase 


Sink 


Characteristics 


Paramagnetic I 




A" ^ 

A ^ -oo 


Low Concentration 
Nonmagnetic impurities 
Dilute Sublattice II 


Ferromagnetic I 


J +0O 

K ^ -oo 
A -oo 


Mag. Order 


Paramagnetic II 


J^O 
K ^0 

A +00 


High Concentration 
Nonmagnetic impurities 


Ferromagnetic II 


J ^ +0O 

A" ^ -oo 

A +00 


Mag. Order 



TABLE I: Phases and Corresponding Sinks 



V. RESULTS 

The results of our investigation, detailed below, probed 
the effects upon ordering of varying the temperature 
(1/J) and vacancy concentration ( A/ J) for a series 
of planes of constant biquadratic coupling K/J. In 
each plane in parameter space, bulk phases have been 
mapped out by exhaustively analyzing the nature of 
renormalization-group trajectories that arise for various 
initial starting parameters (J, if, A). This exploration 
reveals three qualitatively unique regions (or basins of 
attraction), each sharing renormalization-group trajec- 
tories that flow to common sinks, as detailed in Table 
I. 

The two paramagnetic phases arising. Dense Param- 
agnetic and Dilute Paramagnetic, are distinguished from 
one another via the flow of the crystal-field interaction 
term. A flow to -oo corresponding to a high density of 
occupied sites, whereas a flow to +oo corresponds to a 
system dominated by nonmagnetic spin sites and much 
too dilute for magnetic ordering to occur. A Dense Ferro- 
magnetic, complete with crystal-fleld interaction flowing 
to -oo, also arise in phase space, with a renormalization 
group flow that sees the bilinear (J) interaction flowing 
to -|-oo, and, the biquadratic (K) coupling flowing to -oo 
. The ferromagnetic phase can be reached by decreasing 
the temperature, or, by decreasing the concentration of 
nonmagnetic impurities. 

In the plane of constant biquadrupolar coupling 
K/J — 1 we find each of the phases detailed above. 
The Dilute Paramagnetic phase exists at all temperatures 
for larger, positive values of the crystal-field interaction 
(A/ J); the lack of magnetic order consistent with the ex- 
pectations for a lattice dominated by nonmagnetic impu- 
rities. At high temperatures (1/J > 70) we find the Di- 
lute Paramagnetic phase separated from the Dense Para- 
magnetic phase via a first order phase boundary. This 
phase boundary is entirely vertical thus it is traversed 
only via a changing crystal-field interaction. The first 
order boundary separating the two paramagnetic phases 
terminates at a high temperature critical point C. 
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FIG. 3: Parameter space, with K/J = 1 and 
{p,mi,m2,pA,PB) = (4, 8, 9, 40, 1), depicting different 
basins of attraction and associated phases with critical 
endpoint (E) and critical point (C). Solid lines represent 
second-order transitions, whereas dashed lines represent 
first-order transitions. 
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FIG. 4: Parameter space, with K/J — 5 and 
{p,mi,m2,PA,PB) = (4, 8, 9, 40, 1), depicting different 
basins of attraction and associated phases with critical 
endpoint (E) and critical point (C). Solid lines represent 
second-order transitions, whereas dashed lines represent 
first-order transitions. 
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FIG. 5: Parameter space, with K/J — 10 and 
{p, mi, 1712, p A, Pb) = (4, 8, 9, 40, 1), depicting different basins 
of attraction and associated phases with critical endpoint (E) . 
Solid lines represent second-order transitions, whereas dashed 
lines represent first-order transitions. 



FIG. 6: Parameter space, with K/J — 15 and 
{p,mi,m2,pA,PB) = (4, 8, 9, 40, 1), depicting different basins 
of attraction and associated phases with critical endpoint (E) 
and critical point (C). Solid lines represent second-order tran- 
sitions, whereas dashed lines represent first-order transitions. 



In this same plane, K/J — 1, the Dense Paramagnetic 
phase is driven to magnetically order to a Ferromagnetic 
phase with a decrease in temperature via second order 
transition at approximately 1/J ^ 70. This line of crit- 
icality terminates upon intersection with the line of first 
order transitions, at critical endpoint E. The Ferromag- 
netic phase persists for the more negative values of the 
crystal-field interaction, corresponding to a low concen- 
tration of nonmagnetic impurities on the lattice allowing 
magnetic order to propagate. An increase in the crystal 
field drives the system to disorder to the Dilute Param- 
agnetic phase via a first order transition. This line of 
first order transition curves at low temperatures allow- 
ing for the possibility of the Dilute Paramagnetic phase 
to magnetically order with a decrease in temperature for 
a fixed concentration of occupied sites, e.g. A/J = 0.8. 

An increase of the biquadratic coupling to K/J ~ 5 
results in a shift of the first order line of transitions to 
larger A/J. This first order phase boundary separates 
the two paramagnetic states at high temperatures, and 
the ferromagnetic and dilute paramagnetic phases at in- 
termediate and low temperatures. The low temperarture 
phase boundary no longer has the curved feature allow- 
ing for the ordering of the Dilute Paramagnetic phase 
with a decrease in temperature, as seen in the K/ J — \ 
plane. Transition from the Ferromagnetic phase to the 
Dilute Paramagnetic phase driven only with a change in 
the concentration of the nonmagnetic impurities. 

Further increase of the biquadratic coupling to K/ J ~ 
10 reveals a phase diagram that is qualitatively the same 
complete with critical point (C) and critical endpoint (E) 



topologies. However, at low temperatures the boundary 
separating the Ferromagnetic and Dilute Paramagnetic 
phases has curved in the opposite direction as originally 
observed (in the K/J — 1 plane). This phase diagram 
revealing the possibility, at certain concentrations (e.g. 
A/J = 5.6), that the Ferromagnetic phase can surpris- 
ingly disorder to the Dilute Paramagnetic phase with a 
decrease in temperature. 

In the last two phase diagrams considered, K/J = 15 
and K/ J = 20, additional increase in the clustering bias 
(or biquadratic coupling) drives the line of first order 
transitions to larger crystal- fields (A/J) with an increas- 
ing curvature at low temperatures. Interestingly, in both 
of these planes, it is possible to magnetically order the 
Dense Paramagnetic phase via a decrease in temperature 
into the Ferromagnetic phase. And, with an additional 
decrease in temperature, at certain concentrations (e.g. 
A/J — 11.0 in the K/J — 20 plane), it is possible for the 
Ferromagnetic phase to disorder to the Dilute Paramag- 
netic phase. An interesting reentrant phenomena albeit 
to two different paramagnetic states. 

In each plane (of constant biquadratic coupling) con- 
sidered in the present investigation, the high tempera- 
ture transition from the ferromagnetic phase to the dense 
paramagnetic phase was found to be second-order. This 
critical line has been probed and critical exponents calcu- 
lated by linearizing the recursion relations, as discussed 
in Section 4, while maintaining five scaling fields associ- 
ated with J, K, A, iJ, and L. Associated with the scaling 
fields J, K, and A calculation of the eigenvalues of the re- 
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FIG. 7: Parameter space, with K/J = 20 and 
{p, mi, 1712, PA, Pb) = (4, 8, 9, 40, 1), depicting different basins 
of attraction and associated phases with critical endpoint (E) 
and critical point (C). Sohd hnes represent second-order tran- 
sitions, whereas dashed Unes represent first-order transitions. 

suiting recursion matrix yields: two relevant eigenvalues, 
A2 = 19.6 and A4 — 2.00, corresponding to critical scal- 
ing exponents of 1/2 = 4.29 and 7/4 — 1.00, respectively; 
and, an irrelevant eigenvalues with Ag = —1.63. A simi- 
lar analysis associated with scaling fields H and L yields 
two relevant eigenvalues: Ai = 8.00 and A3 = 2.00, cor- 
responding to critical scaling exponents of y2 = 3.00 and 
2/4 = 1.00. 

VI. SUMMARY 

This paper considers a dilute Ising ferromangnet with 
annealed vacancies and attractive biquadratic interac- 



tions. Phase diagrams have been calculated while varying 
the temperature and concentration of annealed vacancies 
in the system while maintaining constant, attractive bi- 
quadratic couplings. These results have been produced 
using renormalization group analysis with a hierarchical 
lattice. Critical exponents have been calculated and each 
basin of attraction interpreted. 

This study has considered several planes of constant, 
attractive biquadratic couplings. Three basins of attrac- 
tion were found each corresponding to a unique bulk 
phase in parameter space. Two paramagnetic regions 
(Dense Paramagnetic and Dilute Paramagnetic), and a 
Ferromagnetic region. First-order boundaries separate 
the Dense Paramagnetic/Dilute Paramagnetic regions at 
high temperatures; this boundary terminating at a high 
temperature critical point C. The Ferromagnet region is 
separated from the Dilute Paramagnetic region also with 
a first-order line of transitions. The high temperature 
Dense Paramagnetic/Ferromagnetic second-order line of 
transitions terminates at a critical endpoint E. An in- 
crease in the biquadratic coupling results in similar phase 
diagrams with greater curvature introduced in the low 
temperature regime thus allowing transitions between the 
Ferromagnetic and Dilute Paramagnetic regions with a 
change in crystal-field interaction and/or temperature. 
The results presented in this paper will provide further 
insight into the role of nonmagnetic impurities upon or- 
dering transitions in ferromagnetic materials as we vary 
the density and clustering bias in the system. 
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